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Abstract 

We investigate the cooperative effects on optical forces in a system of N 
two level atoms confined to a volume of dimension less than A^, where A is 
radiation wavelength and driven by a coherent radiation field with a spatial 
profile like Laguerre-Gaussian or ideal Bessel beam. We show a dramatic 
enhancement on optical forces as well as the angular momentum imparted 
to the atom by a factor of . 

PACS number(s) 42.50 FX , 32.80 Pj 



Recently, there is a considerable interest in cooling and trapping of alkali 

atoms at low temperatures and subsequent observation of 

Bose-Einstein (BEC) condensation [1,2] which has a property of a gen- 
erating laser like source of matter waves. As demonstrated by S. Inoyu et. 
al, these neutral matter waves confined in a cigar shaped BEC [3] show a 
highly directional super radiant Rayleigh scattering when pumped by an off 
resonant laser bcam[4]. This is attributed to a cooperative effect due to 
N atoms spaced to within a wavelength of the laser and is similar (though 
not identical) to Dicke's superradiance[5] phenomenon. It required nearly 
twenty years to observe the superradiant intensity which is proportional to 

accompanied by a subsequent reduction in life time [6]. There has since 
been a resurgence of Dicke's idea after its observation in atomic system of 
two ions by DeVoe and Brewer [7]. They observed a transition rate of two 
ions (r(i?)) as a function of separation(i?) and found that T{R)/To ~ 2 when 
R ^ 0; where Fq is transition rate of a single isolated ion [7]. Very recently, a 
clear experimental demonstration of optical superradiance is demonstrated 
in a thin material sample by Greiner et. al.[8].s 

Dicke's superradiance is observed in many electronic systems of lower di- 
mension too. The systems studied are quantum dots by Brandos et. al.[9], 
and theoretical possibility of the Dicke superradiance is predicted by the 
same group in semiconductor heterostructure[10] where electrons and holes 
in lowest Landau levels act as two level atom— like system. Dicke's superra- 
diance has been proposed as an effective means to reduce spontaneous decay 
time of an excited atom to cool and manipulate atoms with narrow width 
by Djotyan et. oL[ll]. 

In many applications of condensate matter waves, one needs to handle 
the motion of a single atom [12] or of an atomic built— up in a trap[21] or 
amplification of input matter wave[13]. Individual atomic motion can also 
be guided either in free space or in a metallic or dielectric interface[14] or in 
the wave guidc[15]. Most of these studies are confined to manipulate single 
atom in external coherent field. Of late, Gangl and Kitsch have theoretically 
shown that an ensemble of two level atom, coupled to high Q cavity mode 
could be cooled by considering the interaction between the cavity field and 
the motion of atoms in the ensemble[16] indicating an interest in modifica- 
tion of properties due to collective effect. In the light of the resurgence of 
interest in studying collective effects like superradiance, it would be inter- 
esting to investigate the nature of optical forces in a collective system of N 
atoms in Dicke state interacting with a common radiation field. The aim of 
this letter is to show that cooperative effects modify optical forces radically. 

As a model of cooperative effects we consider a collection of N identical 
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two level atoms of transition frequency ujo confined to a volume less than 
where A is a wavelength of transition and driven by a coherent c.w. 
laser field of frequency ujl. Such a model was first considered by Dicke[5] to 
explain superradiance. In a typical superradiance phenomenon, the intensity 
of emitted radiation by the N atom system scales as N^. This implies that 
a collective decay coefficient of where 27 is Einstein's A coefficient. We 
show in this letter that the optical forces (dipolar and dissipative) on the 
above A*" atom system also scale as thereby justifying the title of this 
letter. We also show that the maximum angular momentum imparted to 
the system by the beam is N'^lh where / is quantum number associated with 
the coherent field. 

The equation of motion for the reduced atomic density operator p for 
N identical two level atoms confined to a single site and driven by a single 
mode c-w laser of frequency lol is, 

dp 

ih— = [H, p] + in'j{2'Ki2p7r2i - P7r2i7ri2 - 7r2i7ri2p) (1) 

where, 

for ct*^ atom with i and j taking values 1 and 2 for level | 1 > (ground state) 
and I 2 > (excited state). 27 is Einstein's A coefficient. The total tt operator 
is TTjj = Y^a=i ^ij^ has an algebra 

[7r21,7ri2] = 7r22 - TTii = 27r3 
^3, <2] = -"<2 

[7r2i,7rr2] = -m{m - 1X2"' + 2m7riri7r3 

[7ri2,7r^i] = -n(n - l)7r^f ^ - 2n7r^f 
Hamiltonian of atoms and of field is 

H = tiu;oTr22 - ih{7T2iaf{r)e'^-'^^^ - h.c.) (2) 

It is convenient to make a transformation to the frame of reference of 
the applied coherent field by defining, 
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p = exp(ia;it7r22) pexp(-za;it7r22)- 
The equation(l) is then changed to, 



^ = -^A[7^22,p]-[7^2l«/(r)e*'^-7^l2a7*(r)e-^^p] 

+7(2vri2pvr2i - 7r2i7ri2p - pvr2i7ri2) (3) 

The steady state solution pgs of equation(3) is defined by ^ = and 
has a form [17] 

1 ^ 

/'- = :^ E «mn(-5)-'"(-/)-"vrM (4) 

m,n=0 



where, 



r{m + 5+l)T{n - -5 + 1) 

r(n + i)r(m + i)r(i + s)r{i - s) ' 



(5) 



and 



^ ^ (jV + m + l)!(m!)2 



A = a;o-a;L, 5 = — , g 



7 7 

Calculations of average values of dipolar as well as dissipative forces 
involve averages {7^21) and {7^12) where, 

< TTij >= Trip-Kij). (7) 

The trace is taken over Dicke states 

where k lies between -N/2 and -l-N/2 and k is appropriately summed. The 
expressions for dipolar force < Fdipoie > and the dissipative force < Fdiss > 
are, 

< ^dipole >= ihV I a/(r \){e~'^ < 1,21 > -e^^ < ttu >) (8) 

and 
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< Fdiss >= M a/(r) \Ve{e-'^ < ttsi > +e'^ < tt^ >) (9) 

After evaluating the requisite averages in equations (8) and (9), we obtain 
expressions for forces as, 



{^dipole) 



2hA (V I a/(r)) | 
7 D 

mT{m + S)T(m,-S){N + m+ 1)1 ( \ af{r) ly^^""'^^ 
(2m + l)!(iV - m)!r(l - 5)T{1 + 5) [ j J 



2n I a/(r) I yg 
D 

^ m^r{m + 5)r{m-S){N + m + iy. f | a/(r) , , , 



-(2m-l) 



and 



r(l - S)T{1 + S){2m + 1)\{N - m)! I 7 

It is easy to see that the expressions (10) and (11) reduce to the forces on a 
single atom in an e.m. field when we take N=l as[12] 

< Fa^,ole >= 2n I a/(r) | (V | a/(r) l) ^2 + 2|a/(r) P W ^^^^ 



and 



,2 vg 



< F,.., >= 2^7 I am 1^ ^.^,|^^(,) p^^. (13) 



We have plotted in Fig.(l) 

r (Fdipole) 

and 

n (Fdiss) 

/diss 



[2^752Vg] 

as a function of N for parameters \ S \= 1, g = 1.5 and observe that fdiss 
and fdip tend to a constant value for large N. It means that the sums in 
equations (10) and (11) and the denominator D behave similarly for large 
N. 

Also fdip being a difference of two terms(see eq.8) shows a peak whose po- 
sition is almost insensitive to N for given parameters but its value decreases 
with increasing detuning 6. Note that fdiss is a dimensionless quantity which 
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tends to a value 1 ,an asymptotic limit.We have plotted similiar graphs for 
different values of the parameters | S \ and g and observe that the smaller 
the value of the Rabi frequency, smaller is a value of N around which fdiss 
tends to one. It is natural to scale the quantities like f^ip and fdiss by 
rather than ^ since the width of super radiant pulse is N7 and the peak 
intensity behaves as N'^j [18]. We then see that < Fdipoie > and < F^iss > 
behave as N'^ in the reduced variables g/N and 

For large N, asymptotic analysis of equations (6), (10) and (11) can be 
made following the methods developed by Lawande et al [17]. These authors 
discuss asymptotic expansion of general operator averages of the type 
< 7r2i^(7r22 — 7rii)''7ri2'^ >. Operator averages involved in our case are much 
simpler, namely, < 1121 > and < 7ri2 >. The method, even though exact 
is complicated. The results obtained thus can be obtained in much more 
simple fashion by mean field method. We write 7ri2 = 7r_, 7r2i = 7r+ and 

= (7r22 — 7rii)/2 and write equation of motion for < tt-i- > and tt^. The 
equations involve averages like < 7r±7r2 > which we approximate in mean 
field as, 

< TTiTTz >~< 7r± >< TTz > (14) 

We then go over to reduced variables m± = ruo = ^ff^, t = 

N^t, jJL = ^"j^^'"^^ and = The moment equations in the mean field 
approximation are, 

— — = — — (m+ + m_ j — 2m_)_m_ 
dr 2 

dm+ ium-L. 

— — = /xmo h 2m+mo 

UT I 

dm- ivm^ „ 

— — = iivrio H h 2m-mo (15) 

ar 2 

Using steady state solutions of equations (15) and equations (8) and (9), we 
get the dipole and dissipative force per particle to be 

< Fdipoie > _ n{N-f)u^i{Vii) 
N Ari{l + e) 

and 

N 2(1 + ^2) ^ ^ 

where, = ^[(VFn^) + /3] and r?^ = i[(V/3^ + zy^) - (3] and 4/3 = 
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We now consider an angular momentum imparted to the system of atoms 
by the beams that carry angular momentum such as Laguerre-Gaussian 
(LG) beam or an ideal Bessel beam. The torque imparted will be due to the 
dissipative force. With as a Guoy phase for the beams, we see[19], 



ve 



krz 



+ 



2 
R 

kr'^ 



I 

Cr + -i 



■(1 



2{z^ + zl 
for LG beam and, 



= [^]e^ + {k- §)e, 
for ideal Bessel beam. The torque imparted is. 



2^2 ^^^^{2p + l + l)zR 



z^ + z% 



z^ + zl 



(18) 



(19) 



|< r >| = |< r X Fdiss >| 
The torque imparted by either beams is. 



< T > 



(20) 



(21) 



2(1+ a 

Equations (16), (17) and (21) arc principal results of this letter. It is clear 
that there is a cooperative effect on forces as well as on torque of an exter- 
nal beam if the number of atoms N are within some coherence distance like 
A. Then such a system of atoms experiences a large force and a very large 
torque and very large angular momentum-justifying the part title Dicke's 
bullet. 



Possibility of observing Dicke's bullet is best in gaseous systems.Even 
though, the systems of transition metal ions in solid state matrix as used 
by DeVoe and Brewer [7] and by Greiner et al [8] show superradiance in 
most clear terms it probably is difficult to observe Dicke's bullet in these 
systems . In practice , it is impossible to arrange N atoms at a point. 
However, the correlations between atoms within a distance '^'^^ "^^^^ 
give rise to superradiance. Gaseous samples of sufficiently low density ensure 
negligible interaction between the dipoles . Among many gaseous systems, 
superradiance in Tl-Hg [20] and Cs atomic beam was studied sometime 
ago. If a Laguerre -Gaussian (LG) beam of superradiant frequency falls 
on such a system in a direction opposite to the existing pulse immediately 
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after extinction of the exciting pulse , then, the superradiant atom will 
experience a force as given by Eq. (17) and subsequent acceleration a .The 
cm. velocity imparted to the radiating atoms is .This will Doppler 

shift the peak frequency of u of superradiant transition to {u — which 
may be observable. 

In conclusion , we have shown that the dissipative force and the torque 
for the Dicke system scale as N'^ if N atoms , radiating cooperatively, arc 
exposed to a coherent radiation field with spatial profile such as LG beam 
or ideal Bessel beam. 
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Figure Caption 

Figure 1 shows fdip and fdiss as a function of number of atoms N. ff^p 
and fdiss are dimensionless as defined in the text. 
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Fig. 1 PVP and SVL 
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